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ABSTRACT 

We present results from 30 nights of observations of the open cluster NGC 7789 with the 
WFC camera on the INT telescope in La Palma. From ~900 epochs, we obtained lightcurves 
and Sloan r' — i' colours for ~33()00 stars, with ~2400 stars with better than 1% precision. 
We expected to detect ~2 transiting hot Jupiter planets if 1% of stars host such a companion 
and that a typical hot Jupiter radius is ~ 1.2 Ilj. We find 24 transit candidates, 14 of which we 
can assign a period. We rule out the transiting planet model for 21 of these candidates using 
various robust arguments. For 2 candidates we are unable to decide on their nature, although 
it seems most likely that they are eclipsing binaries as well. We have one candidate exhibiting 
a single eclipse for which we derive a radius of 1.8 Ao'qo-Rj- Three candidates remain that 
require follow-up observations in order to determine their nature. 

Key words: planetary systems - methods: data analysis - Hertzsprung-Russell (HR) diagram 
- binaries: eclipsing - open clusters and associations: general - open clusters and associations: 
individual: NGC 7789 


1 INTRODUCTION 

The surprising existence of short period (~4 days) Jupiter mass 
extra-solar planets (termed “hot Jupiters”) confirmed by radial ve¬ 
locity measurements in the last decade has shown us that plane¬ 
tary systems exist in patterns unlike that of our own Solar System. 
The class of hot Jupiter planets (P < 10 days and M sin i < lOAfj) 
makes up ~19% (27 out of 140 as of 01/02/05) of the planets 
discovered to date (Schneider 1996) and ~1% of ne arby solar 
type stars host such a companion iButleret h] hooch Recently 
we are starting to see the fruits of current transit surveys. OGLE 
have produced over 100 transit candidates du ring two sea sons 

llJdalski et all2002aHUdalski et all2002bl:IUdalski et all2003ll . by 

far the most prolific transit survey. EXPLORE/OC have produced 
a handful of transit candidates that are currently being followed 
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up spectroscopically iMallen-Omelas et aTll2003l) and a search of 
the MACHO phot ometry database has revealed nine transit can¬ 
didates iDrake & Cookll20()4ll . To date there a re seven confi rmed 
transi t ing extra-solar pla nets, HD 209458b l Charb onneau et all 
[2000: Brown et al l 1 2001 ) discovered first by the radial veloc¬ 
ity method, and OGLE-TR-56bjKonadd eydj2003|). OGLE-TR- 
113b iBouchvet al|2004l: iKonacki et al Jl2()04T>. OGLE- TR- 132b 
Bouchv et al JI2004I) . OGLE-TR-11 lb Pont et alJl2004. TrE S-1 
Aions^TaiTzOCMl) and OGLE-TR-10b Konack^^T2004bl) dis¬ 
covered first by the transit method. The spectroscopic follow-up of 
the OGLE transit candidates has revealed a new class of short pe¬ 
riod planets called the “very hot Jupiters”. Such planets have pe¬ 
riods less than the 3 day cut-off identified in the sample of radial 
velocity planets. 

In order to discover transiting hot Jupiters through a photomet¬ 
ric survey, one requires high cadence, high accuracy observations 
(< 1% accuracy per data point with a duty cycle > 8 data points 
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per hour) of many stars (> 10 4 ) simultaneously over long observ¬ 
ing runs (> 10 nights). Any transit candidates (stars that show at 
least one eclipse event) may be subsequently followed up by ra¬ 
dial velocity (RV) measurements in order to determine the com¬ 
panion mass, or at least an upper limit to the mass. However, due 
to the long integration times on large telescopes required for RV 
follow-up and the high frequency of planetary transit mimics, it is 
prudent to try and rule out the transiting planet model for as many 
transit candidates as possible via simple supplementary observa- 
tions and/or further analysis of the lightcurve l Charbonneau 2003; 


ISirko & Paczvhskil2003i:lDrakjl20()^ICharbotmeail2004fr 

The study of open clusters for transiting planets has a num¬ 
ber of advantages over fields in other parts of the sky or Galactic 
plane. While providing a relatively large concentration of stars on 
the sky (but not so large as to cause blending problems as in the 
case of globular clusters observed from the ground), they also pro¬ 
vide a set of common stellar parameters for the cluster members. 
These are metallicity, age, stellar crowding and radiation density. 
Also, the fainter cluster members are smaller stars and therefore 
they are likely to show deeper transit signatures, helping to off¬ 
set sky noise contributions. The identification of the cluster main 
sequence in the colour magnitude diagram allows the assignment 
of a model-dependent mass and radius to each photometric cluster 
member, and assuming a law relating extinction to distance for the 
field allows the assignment of a model-dependent mass, radius and 
distance to all stars in the field under the assumption that they are 
main sequence stars. Transit candidates with well defined phased 
lightcurves may therefore be analysed in detail as to whether they 
are consistent with a transiting planet model. An estimate of the 
fraction of stars hosting a hot Jupiter (referred to as the hot Jupiter 
fraction) may be obtained by comparing the number of hot Jupiters 
that are actually detected to how many one would expect to detect 
using the knowledge of the star properties and the lightcurves them¬ 
selves. The dependence of the hot Jupiter fraction on the cluster 
parameters may then be investigated by extending the experiment 
to other open clusters. 


transiting planets in open clusters l von Braun et al 

2004a). 

These include UStAPS (Street et al. 2003; this paper; 

;H 

ood et al. 

2005 

). EXPLORE/OC iMallen-Omelas et alJ 2003 


Lee et al. 

2004 

; von Braun et al. 2004b), STEPSS 1 Burke etal 

2004) and 

PISCES (Mocheiskaet al. 2002; Mocheiskaet al. 2004 1 . A num- 


ber of transit candidates have been put forward by these groups but 
none have been confirmed as transiting planets so far. In this sur¬ 
vey, as with open cluster transit surveys in general, we are photo¬ 
metrically observing faint stars (16 to 21 mag). This makes follow¬ 
up observations difficult but not impossible. However, many of the 
stars observed are of later spectral types than those probed by the 
RV surveys since RV surveys are limited to bright solar neighbour¬ 
hood FGK stars. In particular, our survey of NGC 7789 probes the 
spectral types F5 to M5 (see Fig. 2) including a large proportion 
of K and M stars. Furthermore, we are searching for planets out to 
distances well beyond the solar neighbourhood. 

The observations of open cluster NGC 7789 is the subject of 
this paper. The main parameters of the cluster are shown in Table 
1. For a_good review of previous relevant work on this cluster see 
iGim et aL il l998 1. In Section 2 we report on the observations made, 
in Section 3 we present in detail the data reduction/photometry, in 
Section 4 we present the astrometry and colour data, in Section 
5 we describe the transit detection algorithm and in Section 6 we 
present a detailed analysis of the transit candidates. Finally, in Sec¬ 
tion 7 we outline our conclusions and future work. 


Table 1. Properties of the open cluster NGC 7789. Data taken from 
http://obswww.unige.ch/webda by Mermilliod, J.C. and the SIMBAD 
database. 


RA (J2000.0) 

23 h 57 m 

Dec (J2000.0) 

+56°43' 

l 

115?48 

b 

-5? 37 

Distance (pc) 

2337 

Radius 

~16 / 

Age (Gyr) 

1.7 

[Fe/H] 

-0.08 

E(B - V) 

0.217 


2 OBSERVATIONS 

We observed the open cluster NGC 7789 (see Table 1) using the 
2.5m Isaac Newton Telescope (INT) of the Observatorio del Roque 
de los Muchachos, La Palma, in the Canary Islands during three 
bright runs with dates 1999 June 22-30, 1999 July 22-31 and 2000 
September 10-20. For brevity, these runs shall be refered to from 
now on as 1999-06, 1999-07 and 2000-09 respectively. We used 
the Wide Field Camera (WFC) which consis ts of a 4 EEV CCD 
mosaic where each CCD is 2048x4096 pixels IWalton et all200lh . 
The pixel scale is 0.33”/pix and field of view ~0.5°x 0.5°. The 
gain and readout noise values for each chip were calculated au¬ 
tomatically during the preprocessing stage of the data reduction 
(see Section 3.1). The mosaic field was centred on NGC 7789 at 
a = 23 h 57 m 30 s and S = +56°43'41". 

The usual procedure for each night was to obtain ~5 bias 
frames and ~8 sky flat frames at both the beginning and end of the 
night. Observations on NGC 7789 in the runs 1999-06 and 1999-07 
consisted of ~6 pairs of 300s exposures taken every ~50 minutes 
during the later part of each night. Observations in the 2000-09 run 
consisted of sequences of ten consecutive 300s exposures followed 
by a bias frame, repeated continuously throughout the whole of 
each night. With a readout time of 100s and various losses due to 
bad weather/seeing and telescope jumps, this resulted in a total of 
880 x 300 s exposures in Sloan r' over the three runs, with 691 of 
these exposures from the 2000-09 run alone. During the 2000-09 
run, we also took 5 images of NGC 7789 with varying exposure 
times in Sloan i !, along with 5 sky flat frames, in order to provide 
us with the necessary colour information. 


3 CCD REDUCTIONS 

3.1 Preprocessing: CCD Calibrations 

Each run and each chip was treated independently for the purpose 
of the reductions. The reduction process was carried out by a single 
C-shell/IRAF script that runs according to a user-defined parame¬ 
ter file. Bad pixels were flagged in a user-defined detector bad pixel 
mask, and ignored where relevant. The reductions consisted of the 
standard bias subtraction and flat fielding of the science frames us¬ 
ing appropriate master calibration frames. The gain and readout 
noise of each chip were also determined during the reductions. 


3.2 Photometry: Difference Image Analysis 

Differential photometry on the reduced science frames in the Sloan 
r' filter was accomplished using the method of difference im- 
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INT/WFC NGC7789_INT_Chip1r 8631 Stars at 607 Epochs 


INT/WFC NGC7789_INT_Chip2r 7625 Stars at 604 Epochs 




star magnitude 


star magnitude 


(a) CCD 1 


(b) CCD 2 


NT/WFC NGC7789JNT_Chip3r 8411 Stars at 591 Epochs 


NT/WFC NGC7789_INT_Chip4r 8830 Stars at 613 Epochs 



star magnitude 
(c) CCD 3 



Figure 1. Plots of standard deviation (RMS) of the lightcurves against mean instrumental Sloan r' magnitude for all stars from each chip for the 2000-09 run. 
The lower curve in each diagram represents the theoretical noise limit for photon and readout noise. 


age analysis iAlard & LuDtoiJfr99&lAlardil200dh . Our implemen¬ 
tation of this procedure was_adapted from the code written for the 
MOA project iBond et alJ200lt) . and it consists of three automated 
scripts. The process is described below in relative detail since the 
same procedure has subsequently been used on many other data sets 
(for example I Arellano Ferro et all2004 . Bad pixels are ignored in 
the operations that the scripts perform. 

(i) The first script constructs a reference frame from selected 
frames with good seeing, and a star list from the reference frame. 
First, stars are detected and matched between the best seeing frames 
in order to derive a set of linear transformations and geometrically 
align the frames. The frames are then combined into a mean ref¬ 
erence frame using the exposure times of the individual images as 
weights. The reference frame is analysed using IRAF's DAOPhot 
package jStetsonllll987t) . The package identifies stars on the refer¬ 
ence frame and chooses a set of 175 point-spread function (PSF) 
stars. A “penny2” PSF function that varies quadratically with po¬ 
sition, along with a lookup table of residuals, is solved for. The 
neighbours of the PSF stars are then subtracted using this solution, 


and a new PSF function is solved for. This new solution is used 
to measure the instrumental fluxes and positions of all stars on the 
reference frame. The result is a reference frame with a correspond¬ 
ing star list. We used 13 consecutive best seeing images (~l") to 
construct the reference frame. 

(ii) The reference frame is used to produce a set of difference 
images. The main idea behind difference image analysis is that an 
image frame I(x,y) is related to the reference frame R{x,y ) via 
the following equation: 

I(x,y) = R® K(x,y) + B(x,y) (1) 

where 

R® K(x,y) = J J R(x — u, y — v) K(u, v, x, y) du du (2) 

Here B{x,y ) represents the change in the sky background and 
K(u, v, x, y) is a spatially-varying convolution kernel relating the 
point-spread function on the reference frame to the point-spread 
function on the image frame at spatial position x, y. We model the 
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Table 2. Magnitude offsets Ami and Am .2 added to the lightcurve data 
points from the 1999-06 and 1999-07 runs respectively. 


Chip No. 

Ami 

&1 

A m 2 

<7 2 

1 

-0.832 

0.037 

-0.753 

0.041 

2 

-0.698 

0.036 

-0.843 

0.043 

3 

-0.596 

0.031 

-0.404 

0.054 

4 

-0.527 

0.036 

-0.470 

0.037 

Run: 

1999-06 


1999-07 



convolution kernel: 

K(u,v,x,y) = Oi(x,y)bi(u,v) (3) 

i 

as the sum of a set of basis functions bi(u,i>) each formed as a 
product of a two-dimensional Gaussian function of u and v with a 
polynomial of degree 2 in u and v. For the basis functions we use 
3 Gaussian components with sigmas of 2.1 pix, 1.3 pix and 0.7 pix 
and associated polynomial degrees of 2, 4 and 6 respectively. To al¬ 
low for the kernel’s spatial dependence, the coefficients di(x, y ) are 
polynomials of degree 2 in * and y. The kernel is also normalised 
to a constant integral over u and v for each x and y, thus ensuring a 
constant photometric scale factor between the reference frame and 
image frame. We model the differential sky background B(x, y) as 
a polynomial of degree 2 in x and y. We solve for K(u, v, x, y) 
and B(x,y) in the least-squares sense for each science frame by 
fitting to pixel boxes around selected bright stars distributed uni¬ 
formly across the reference frame. The kernel is assumed to be in¬ 
dependent of x and y within each box. A difference image is then 
constructed for each science frame by rearranging Eqn.Qto the fol¬ 
lowing form and using the solutions for K ( u, v, x, y) and B(x, y): 

D(x, y) = I(x,y) - R®K{x,y) - B(x,y) (4) 

The difference image D(x, y) should simply be an image repre¬ 
sentative of the Poisson noise in I(x,y). However, any objects 
that have varied in brightness in comparison to the reference frame 
should show up as positive or negative pixel areas on the difference 
image which may be measured to obtain the differential flux. In 
our analysis, each chip was split up into 8 square sections and the 
difference image constructed from solving for the kernel and differ¬ 
ential sky background in each section. Also, a high signal-to-noise 
empirical PSF for the reference frame is constructed in each section 
by stacking up a set of stamps centred on suitable bright stars. 

(iii) The third script measures the differential flux on each dif¬ 
ference image via optimal PSF scaling at the position of each star. 
The normalised and sky-subtracted empirical PSF constructed for 
each square section of the reference frame in step (ii) is convolved 
with the kernel corresponding to the current difference image sec¬ 
tion. The convolved PSF is optimally scaled, at the position of the 
current star, to the difference image. A 3<r clip on the residuals of 
the scaling is performed, and one pixel rejected. The scaling and re¬ 
jection is repeated until no more pixels are rejected. The differential 
flux is measured as the integral of this scaled PSF. 

A lightcurve for each star was constructed by the addition of 
the differential fluxes to the star fluxes as measured on the reference 
frame. The following equations were used: 

/tot(*)=/ref+^^ (5) 


m(t) = 25.0- 2.5 log(/ tot (f)) (6) 

where /tot(f) is the star flux (ADU/s) at time f, / re f is the star 
flux (ADU/s) as measured on the reference frame, /djff(i) is the 
differential flux (ADU/s) at time t as measured on the difference 
image, p(t) is the photometric scale factor (the integral of the kernel 
solution over u and v) at time t and m(t ) is the magnitude of the 
star at time t. Uncertainties are propagated in the correct analytical 
fashion. 

Flux measurements were rejected for a x 2 pi x_1 > 5.0 for 
the PSF scaling, and for PSFs with a FWHM>7.0 pix, in order to 
remove bad measurements. Hence, all the stars have differing num¬ 
bers of photometric measurements. In each run, lightcurves with 
less than half of the total possible epochs were rejected. For the 
2000-09 run this analysis produced 8631 lightcurves on Chip 1, 
7625 lightcurves on Chip 2, 8411 lightcurves on Chip 3 and 8830 
lightcurves on Chip 4 (centred on the cluster). Fig. Q shows a di¬ 
agram of the RMS scatter in the lightcurves against instrumental 
magnitude for the 2000-09 run for each chip. Similar diagrams 
were produced for the 1999-06 and 1999-07 runs but are not shown 
here for brevity. 

Since each run was treated independently for the reductions, 
each chip has three different reference frames and hence each star 
has three different reference magnitudes. For a particular star, let 
us denote the reference magnitude from the 2000-09 run minus the 
reference magnitude from the 1999-06 run by Ami and the refer¬ 
ence magnitude from the 2000-09 run minus the reference magni¬ 
tude from the 1999-07 run by Am. 2 . For each chip, we have calcu¬ 
lated the unweighted mean of Ami and Am 2 over all stars on that 
chip. We then added the resulting Ami and Arri 2 to the lightcurve 
data points in the 1999-06 and 1999-07 runs respectively. The val¬ 
ues of the means Ami and Am 2 for each chip along with the stan¬ 
dard deviations about the means <n and 02 respectively are pre¬ 
sented in Table 2. 

As can be seen from Fig.Q we have obtained high precision 
photometry with an RMS accuracy of ~3-5mmag at the bright end. 
Most stars are limited by sky noise because all three runs were dur¬ 
ing bright time. However, the "backbone” of points on each dia¬ 
gram lies above the theoretical limit by a factor of ~1.5-2.0 de¬ 
pending on the chip being considered. We put this down to system¬ 
atic errors in the data due to a subset of low quality difference im¬ 
ages and/or sections of difference images that were produced from 
science frames taken on nights of poor quality seeing/atmospheric 
conditions. 


4 ASTROMETRY AND COLOUR DATA 
4.1 Astrometry 

Astrometry was undertaken by matching 358 stars from the four 
referenc e frames (o ne for each chip) with the USNO-Bl.O star cat¬ 
alogue {Monet e t al|[2003|t using a field overlay in the image dis¬ 
play tool GAIA lDranei&2OQ0l) . The WFC suffers from pincushion 
distortion, hence it was necessary to fit a 9 parameter astromet¬ 
ric solution to the reference frames in order to obtain sufficiently 
accurate celestial coordinates for all the stars. The 9 parameters 
are made up of 6 parameters to define the linear transformation 
between pixel coordinates and celestial coordinates, 2 parameters 
to define the plate centre and 1 parameter to define the radial dis¬ 
tortion coefficient. The starlink package ASTROM 4Wallaceh998l) 
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Table 3. Number of stars with a Sloan r' lightcurve, and the number of such 
stars with a Sloan r' — i 1 colour index. 


Chip 

No. 

No. Stars With 

A Lightcurve 

No. Stars With A 
Lightcurve And An 
r' — i' Colour Index 

Percentage 

1 

8631 

8497 

98.4% 

2 

7625 

7576 

99.4% 

3 

8411 

8290 

98.6% 

4 

8830 

8672 

98.2% 

Total: 

33497 

33035 

98.6% 


was used to do the fit and the achieved accuracy was ~0.4 arcsec 
RMS radially for the 358 matching stars. The astrometric fit was 
then used to calculate the J2000.0 celestial coordinates for all stars 
with a lightcurve. 

4.2 Colour Indices 

The best image in the Sloan i' filter was aligned with the Sloan 
r' reference frame for each chip and the magnitudes of the stars 
were measured using DAOPhot PSF fitting in the same way as they 
were measured on the reference frame in Section 3.2. Table 3 shows 
the number of stars with lightcurves that have Sloan r' — i r colour 
indices as a result. 

4.3 Colour Magnitude Diagrams 

Figs. 2a-d show an instrumental colour magnitude diagram (CMD) 
for each chip. The cluster main sequence is clearly visible. Chip 4 
is centred on the cluster and as expected shows the strongest clus¬ 
ter main sequence. A theoretical cluster main sequence is plotted 
on each diagram over the cluster main sequence. We have used the 
theoretical models of lBaraffe et ail <1998li for the stellar mass range 
0.60 Mq < M* < 1.40Mq, the age of the cluster (1.7Gyr) and 
solar-type metallicity [M/Fl] = 0 in order to predict the main se¬ 
quence absolute magnitudes, colours and radii. Below a mass of 
0.60 Mq the Baraffe model predicts R — I colours substantially 
bluer than the obser ved cluster main sequence, a limitatio n note d in 
iBaraffe et alj|l998 :). As a result we used data from Lang ( 1992 ) for 
the stellar mass range 0.08 Mq < M » < 0.60Mq. The combined 
model for the cluster main sequence supplies an absolute magni¬ 
tude Mr, an absolute magnitude Mi and a stellar radius R * for the 
stellar mass range O.OSMq < M, < 1.40Mg. We interpolated 
this combined model with cubic splines. 

The interstellar medium (ISM) in the Milky Way is mostly 
concentrated in the Galactic plane and the density law governing its 
mean distribution (ignoring small scale variations) can be modelled 
by an Einasto law: 

p(R,z) = p oexp (- exp (-g) (V) 

where R is the Galactocentric distance, 2 is the height above the 
Galactic plane, po is the local density of the ISM, Rq is the dis¬ 
tance of the Sun from the Galactic centre, Hr is the ISM density 
scale height in the R direction and h z is the ISM density scale 
height in the « direction. One may derive the density p of the ISM 
as a function of distance d from the Sun in the direction of the open 
cluster NGC 7789 by using trigonometrical arguments to rewrite R 
and a as functions of d. In this derivation, we have assumed that 


the Sun has Galactic coordinates (R, z) = (8.5kpc, 0.015kpc) and 
th at pq = 0.0 21Mg pc~ 3 , tiR = 4.5kpc, h z = 0.14kpc as given 
in lRobin et alJ (20031). In any wave band, the total extinction A as a 
function of d is proportional to the integral of p(d) over d. Hence, 
absorbing the constant po into a new constant K we have: 

A(d) = K [ p{u) du (8) 

Jo 

Adopting E(B — V) = 0.217 for the cluster (Table 1), we 
calculate the corresponding extinction to be Ar ~ 0.547 and 
Ai « 0.429 in the R and I bands respectively, evaluated with 
a synthetic photom etry co de (XCAL) using a Galactic extinction 
curve from ISeaton t !979 ). This extinction applies to stars at the 
cluster distance d c = 2337pc, and hence, by numerically evalu¬ 
ating the integral in Eqn. [8] we may calculate values for I\ that 
apply to the R and / bands as Kr = 2.20 x 10 _2 magMQ 1 pc 2 
and Ki = 1.73 x 10 _2 magMQ 1 pc 2 respectively. 

We have used the law relating extinction to distance as given 
in Eqn. [8] to correct the absolute magnitudes Mr and Mi of the 
theoretical main sequence to the observed magnitudes R(d) and 
1(d) respectively. In the following equations, the distance d has 
units of parsecs (pc): 

R(d) = Mr + 5 log(d) — 5 + A R (d) (9) 

1(d) = Mi + 51og(d) — 5 + Ai(d) (10) 

where Ar(6) and Ai(d) are versions of Eqn.[8]with K = Kr and 
I\ = Ki respectively. 

Conversions between the Johnson-Cousins R and / magni¬ 
tudes and the Sloan r' and %' magnitudes were done using the fol¬ 
lowing predetermined relations presented on the Cambridge Astro¬ 
nomical Survey Unit (CASU) webpage 1 : 

r = R + 0.275(77 - 7) 4- 0.008 (11) 

r — i! = 1.052(7? — 7) + 0.004 (12) 

Due to the lack of observations of standard stars, it was neces¬ 
sary to fit the interpolated theoretical main sequence to the cluster 
main sequence on the CMD for each chip by eye, after correcting 
for the cluster distance and extinction, via simple r' and r' — i' off¬ 
sets. These offsets are displayed in Table 4. Note that the required 
horizontal and vertical shifts are correlated, since shifts parallel to 
the main sequence would have no effect if the main sequence were 
a straight line. Fortunately, the kink (change of slope) in the main 
sequence near the spectral type K0 (0.8 Mq) allows us to estimate 
both vertical and horizontal shifts. This feature is clearly visible on 
all 4 chips. 

4.4 Stellar Masses, Radii And Distances 

The identification of the cluster main sequence on each CMD al¬ 
lows a model-dependent mass, radius and distance for each star to 
be determined using the theoretical main sequence, assuming that 
each star is a main sequence star. Giant stars (MK Luminosity Class 
III) have absolute magnitudes in the range 1.7 < My ~ Mr < 
—6.5 iLan dl992h . Assuming that the Sun lies in the Galactic plane 
at a distance of 8.5kpc from the Galactic centre (the IAU value) and 

1 http://www.ast.cam.ac.uk/~wfcsur/index.php 
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Instrumental r’-i’ 



Instrumental r’-i’ 


(a) CCD 1 


(b) CCD 2 



Instrumental r’-i' 



(c) CCD 3 


(d) CCD 4 



Figure 2. (a)-(d): Instrumental CMDs for all stars from each chip for the 2000-09 run. The main sequence is visible on each chip, and the theoretical cluster 
main sequence is overlayed as the dashed line. The straight dotted lines are the faint limits for giant stars assuming no extinction and the law relating extinction 
to distance in Eqn. I8l The transit candidates of Section 6 are marked on as solid circles. The errorbars on the right hand side of each diagram represent the 
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Table 4. Offsets determined by eye between observed rf and r' bs — i' b 
magnitudes, and theoretical r'(d c ) and r'(d c ) — i'(d c ) main sequence 
magnitudes for the cluster distance dc. 


Chip No. 

^3- 

1 

Jo 

~~ o 

S- 

Kb 

i 

-0.3 

-0.95 

2 

-0.3 

-0.80 

3 

-0.1 

-0.70 

4 

-0.1 

-0.75 

Estimated 

Error: 

±0.1 

±0.05 


assuming that the Galactic disk has a radius of 14.0kpc i lRohin et alJ 
[1992), then the distance to the edge of the Galaxy in the direc¬ 
tion of NGC 7789 may be calculated as ~8.1kpc using elementary 
trigonometry. The magnitude of the dimmest giant at 8.1kpc as¬ 
suming no extinction is R = 16.2 and, assuming the law relating 
extinction to distance in Eqn.0 the dimmest giant has a magni¬ 
tude of R = 16.9. These faint limits are marked on the CMDs in 
Figs. 2a-d as dotted lines. From this simple argument it can be seen 
that only the brightest stars in our sample will be contaminated with 
giant stars. 

In principle, for each star on the CMD, it is possible to choose 
a value for the distance parameter d in Equations 9 and 10 such 
that the theoretical main sequence passes through the star's posi¬ 
tion on the CMD. The solution d = d* is then the distance to the 
star. The star mass M» and radius R , may subsequently be deter¬ 
mined from where the star lies on the theoretical main sequence at 
a distance of d*. Fig. Et shows the grid of star masses and distances 
used for Chip 4. The solid vertical lines represent lines of constant 
stellar mass (and radius), and are labelled at the top of the diagram 
in units of Mq. The “diagonal” dashed lines represent theoretical 
main sequence models at different distances, and the distances are 
labelled to the right and bottom of the diagram in units of parsecs. 
Fiducial spectral types are marked on the cluster theoretical main 
sequence for clarity. 

Due to the steepness of the theoretical main sequence in the 
CMD for star masses greater than 0.80Mg, the determined star 
properties become more uncertain above O.SOMq. Also, the the¬ 
oretical main sequence that we have used terminates at a mass of 
1.40Mq, which leads to a small region where there are no solu¬ 
tions for d*. In Fig. this region is blueward of the thick con¬ 
tinuous line (corresponding to a mass of 1.40Mq). Stars with no 
solution for d* have masses greater than 1.40Mq (and radii greater 
than 1.70-Rg) and large distances. It is around these stars that it is 
hardest to detect a transiting planet and hence a lack of solution for 
d*, M* and R » will hardly affect the completeness of our survey. 
Table 6 shows the star masses, radii and distances obtained by the 
above procedure for the transit candidates discussed in Section 6. 
The star r’ — i! colours have been corrected where necessary for 
any lightcurve variations (since the reference frame from which the 
r' magnitude was determined has a different epoch to the i frame 
from which the i magnitude was determined). 

In Figs.0i-b we plot the r 1 magnitude and standard deviation 
(RMS) of the lightcurve of each star versus the stellar radius de¬ 
rived from the main-sequence model and the observed r' — i' colour 
index. A vertical stripe of stars is evident with R » ~ 0.75 Rq. This 
arises because of a relatively rapid change in the colour index with 
mass for the theoretical main sequence in this mass range. This ef¬ 
fect is also evident in Fig. 0>, where the vertical iso-mass lines are 


more widely spaced for 0.5Mq < M* < 0.8Mq. If the mass 
function and mass-radius relationship for main sequence stars are 
both smooth, then this effect represents a deficiency in the R — I 
colour index of the stellar models. 


4.5 The Expected Number of Transiting Planets 

In Fig. 03 each point below one of the curves represents a star 
whose lightcurve has sufficient accuracy in our dataset to reveal 
transits by HD 209458b-like planets of the indicated radius. To ob¬ 
tain these detection threshold curves, we note that the signal-to- 
noise ratio for detecting a transit is roughly: 


S/N 


AtovTV 

a 


(13) 


where Am is the transit depth, a is the standard deviation of the 
photometric measurements, and N is the number of data points ac¬ 
quired during the transit. The transit depth in magnitudes is approx¬ 
imately: 


Am 


-2.5 log 



(14) 


where Rp is the radius of the planet and R* is the radius of the 
host star. For random sampling of the orbital period P, the proba¬ 
bility that a given data point catches a transit is At/ P, where At is 
the transit duration. For a HD 209458b-like system (Am « 1.5%, 
At ss 3 h, P ss 3.5 d), this fraction is At/P = 3.5% and thus 
N ~30 of our 880 lightcurve data points would catch a transit. 
For an OGLE-TR-56b-like system (Am S3 1.3%, At S3 1.7 h, 
P S3 1.2 d), At/P = 5.9% and thus N ~52. In Section 5 we 
adopt a conservative transit detection threshold of S/N~10. By re¬ 
arranging Eqn. ll3l and using Ean. 1141 we have plotted in Fig. 03 the 
required RMS accuracy a to detect at S/N=10 a HD 20945 8b-like 
system (N = 30) as a function of stellar radius i?* for planetary 
radii of 0.57?j, 1.0-Rj, 1.4i?j and 2.0i?j. 

In fact, we may count the number of stars in our survey that 
have an RMS accuracy better than that required for a S/N~10 de¬ 
tection with N data points during transit. In Fig. m we plot this 
number against Rp for N = 20, N = 30 and N = 50. The curves 
N = 30 and IV = 50 are representative of a typical hot Jupiter 
(HD 209458b) and a typical very hot Jupiter (OGLE-TR-56b) re¬ 
spectively. In Fig. 01 we plot P(visible) against orbital period where 
P(visible) is the probability that at least M transits have a time of 
mid-transit that occurs during our observations. The values M = 1, 
M = 2 and M = 3 are represented by the upper, middle and lower 
curves respectively. 

We may use Figs. m-d to estimate the expected number of 
transiting planets. Assuming a typical hot Jupiter radius of 1.2i?j, 
then we estimate that there are 3083 and 4625 stars (Fig. 0) whose 
lightcurves have sufficient accuracy in order to detect HD 209458b- 
like and OGLE-TR-56b-like planets respectively. Fig. 01 shows 
that 1 to 3 day planets are almost guaranteed to transit during our 
observations and that 3 to 5 day planets have ~80% probability. 
Assuming that ~1% of stars host a hot Jupiter companion, and that 
~10% of such systems exhibit transits, then we may expect ~2 
stars in our sample to reveal HD 20945 8b-like planetary transits or 
~4 stars in our sample to reveal OGLE-TR-56b-like planetary tran¬ 
sits. In forthcoming work (Bramich & Horne 2005) we report re¬ 
sults of more detailed modelling of the planet detection capabilities 
of our survey based on Monte Carlo simulations that are consistent 
with these results. 
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(C) 



Figure 3. (a): A plot of r' magnitude against star radius for all four chips. The continuous line is the theoretical cluster main sequence, (b): Standard deviation 
(RMS) of the lightcurves against star radius for all four chips. The continuous curves are the detection limits for a HD209458-like planet of the quoted radius 
as a function of star radius (see Section 4.5). (c): Number of stars in our survey that have an RMS accuracy better than that required for a S/NRilO detection 
with N data points during transit as a function of planetary radius, (d): A plot of P( visible) against orbital period for M = 1 (upper curve), M = 2 (middle 
curve) and M = 3 (lower curve). P(visible) is the probability that at least M transits have a time of mid-transit that occurs during our observations. 


5 TRANSIT DETECTION 

We used a matched filter algorithm to search for transits in the 
lightcurves. Adopting a square “boxcar” shape for the transit 
lightcurve, the transit model has 4 parameters: the out-of-eclipse 
magnitude mo, the time of mid-transit to, duration At and depth 
Am. We search for transits with durations ranging from 0.5 h to 
5 h, spanning this range with 12 values of At spaced by factors of 
1.23. We move the transit centroid to through the data in steps of 
At/4. As illustrated in Fig. [4] we fit both a constant and a box¬ 
car transit lightcurve to the data points in a window of width 5Af 
centred on each value of to. Our transit detection statistic is: 


2 _ 2 

q 2 _ Xconst ~ Xtra 
°tra — / 2 

Xout 

Aout - 1 


(15) 


where y/ ra is the chi squared of the boxcar transit fit, Xconst * s ih e 
chi squared of the constant fit, Xout' s ih e c hi squared of the boxcar 
transit fit for the AW out-of-transit data points. The statistic Sp. a 
is effectively the squared signal-to-noise ratio of the fitted transit 
signal renormalised to the reduced chi squared of the out-of-transit 
data points. This modified matched filter algorithm was designed 
to help downweight systematic errors with Xq u1 / (N ou i — 1) > 1 
(and serendipitously, variables), since transit signals should have 
Xout/ (AW — 1) ~ 1- 

The transit detection algorithm outlined above was applied 
to the 1999-07 and 2000-09 runs. Initial tests with ,S\ 2 ra generated 
many spurious transit candidates in which the transit fit matched 
low data points at the beginning or end of a night. To suppress these 
we introduced additional requirements on the number of in-transit 
and out-of-transit lightcurve data points. For the densely sampled 
2000-09 run we required at least 3 in-transit and 8 out-of-transit 
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~a 



Figure 4. An example boxcar transit fit showing the in-transit and out-of- 
transit zones (continuous line) and the constant fit (dashed line). The hori¬ 
zontal axis is time (days) and the vertical axis instrumental Sloan r' magni¬ 
tude. “Statistic” is the value of the transit statistic ,S' 2 , L for this fit. 


lightcurve data points for a transit detection. For the more sparsely 
sampled 1999-07 run we required at least 2 in-transit and 6 out-of- 
transit lightcurve data points for a transit detection. The time sam¬ 
pling in the 1999-06 run was too sparse to support transit hunting 
via the above technique. 

In each lightcurve the highest value of S^ a on each night was 
identified, and those with Sp. a > 100 (equivalent to S/N > 10) 
were retained for closer examination. Table [5] lists for each chip 
the number of raw candidate transits thereby selected over the two 
runs. Despite the high signal-to-noise threshold for detection, 2182 
raw transit candidates were found. A careful visual inspection of 
the corresponding lightcurves lead us to reject the majority of these 
based on a number of criteria. The majority of the raw transit can¬ 
didates (61.8%) were rejected because they appeared to represent 
a single much fainter data point resulting from a “bad” section in 
one of the difference images. Such cases were readily identifiable 
because the lightcurves of many stars triggered a transit detection 
at the same epoch. A large number of variable stars were picked up 
(~ 100 lightcurves = 19.5% of the raw transit candidates), which 
we plan to present in a forthcoming paper. Lightcurves showing 
eclipses with clearly different depths were also assigned as vari¬ 
able stars since a stellar binary is indicated in this case. 

For the remaining raw transit candidates we examined the star 
on the reference frame. This revealed that many of the remaining 
transit signatures were caused by the following (in order of most 
common occurence): 

(i) Image defects detected as “stars” (4.5%). 

(ii) Stars lying on or close to image defects, bad columns and/or 
saturation spikes (4.0%). 

(iii) Stars close to saturation (3.7%). 

(iv) Very closely blended stars (2.1%). 

(v) Stars close to the edge of the CCD (0.8%). 

The reference image for each chip contained a large number of sat¬ 
urated stars along with large saturation spikes which unfortunately 
increased the incidence of such false alarms. 

For the transit candidate lightcurves that survived to this point, 
we checked the difference images for the night(s) of the suspected 
transit(s) by constructing a difference image movie. This revealed 
that a handful of the candidates (0.6%) were the result of a consec¬ 
utive set of poor subtractions at the star position. The other candi- 



10000 


1000 


Reduced Chi Squared 'Out Of Transit’ 


1000 


Figure 5. The transit detection statistic S aa against the out-of-transit re¬ 
duced chi squared Xo ut /(/Voijt — 1). The initial transit detection threshold 
shown by the horizontal line is set at = 100. The strongest transits de¬ 
tected in the lightcurves of the 24 stars that survived subsequent data quality 
tests are plotted as solid circles. The blank semicircular region is saturated 
with test points. 


Table 5. The number of raw transit candidates and remaining transit candi¬ 
date lightcurves after the weeding process for each chip over the two runs. 


Chip No. 

No. Raw Transit 

Candidates 

No. Remaining Transit 
Candidate Lightcurves 

i 

617 

9 

2 

311 

5 

3 

721 

6 

4 

533 

4 

Total: 

2182 

24 


dates clearly showed a flat difference image followed by a growing 
and then diminishing “dimple”, indicating a drop and then recovery 
in the brightness of the star. 

We discuss in Section 6 below the 24 transit candidate 
lightcurves that passed all of the data quality tests outlined above. 
Reference to a transit candidate from now on refers only to one of 
these transit candidate lightcurves. Fig.[5]shows all tests for which 
S t 2 ra >10 and highlights the eclipse with the greatest value of S t 2 a 
for each transit candidate. Table 6 details the number of fully and 
partially observed eclipses that are present for each transit candi¬ 
date and how these eclipses are distributed between the three runs. 
Table 6 also lists the J2000.0 celestial coordinates for each transit 
candidate. 


6 TRANSIT CANDIDATES 
6.1 Theoretical Models 

The lightcurves of the 24 transit candidates selected in Section 5 
were modelled as a star and planet system in the following way. We 
assume spherical stars, a luminous primary of radius R , and a dark 
massless companion of radius Rc in a circular orbit with radius a 
and period P inclined by the inclination i relative to our line of 
sight. The time to is the time of mid-eclipse of the primary by the 
companion. Since we already know 77*, the parameters that need 
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Table 6. Star and lightcurve properties for the transit candidates. The number in brackets indicates the uncertainty on the last decimal place. Columns 2 and 3 
are calibrated r' and r' — i' using the shifts in Table 4. Columns 7-9 indicate the number of fully and partially observed eclipses (f,p) that are present in the 
lightcurve of the corresponding run. The eclipsing stellar binaries are classified into the following categories: E = Eclipsing binary, EA = Algol type eclipsing 
binary, RS = RS Canum Venaticorum type eclipsing binary, CV = Cataclysmic variable. 


Star 

r' 

r' — i ' 

Mass 

Radius 

Distance 

1999-06 

1999-07 

2000-09 

RA 

Dec 

Variable 


(mag) 

(mag) 

(AT©) 

(R©) 

(pc) 

f,p 

f,P 

f,p 

(J2000.0) 

(J2000.0) 

Class 

EB-1 

18.077 

0.569(17) 

0.949(33) 

0.947(42) 

2908(13) 

0,0 

0,1 

0,0 

23 58 36.75 

+56 26 56.1 

E 

EB-2 

19.140 

0.538(9) 

1.028(22) 

1.062(35) 

5762(11) 

0,0 

0,1 

0,2 

23 58 33.88 

+56 37 04.9 

E 

EB-3 

19.346 

0.606(12)“ 

0.897(23) 

0.882(29) 

4346(31) 

0,0 

0,1 

0,1 

23 57 11.92 

+56 31 24.9 

RS 

EB-4 

19.447 

0.999(10) 

0.666(4) 

0.635(4) 

1799(4) 

0,0 

0,0 

0,1 

23 56 46.79 

+56 36 13.8 

E 

EB-5 

20.088 

0.773(18)“ 

0.751(5) 

0.716(5) 

3564(20) 

0,0 

1,1 

1,1 

23 55 58.93 

+56 40 29.6 

E 

EB-6 

17.338 

0.695(45) 

0.759(23) 

0.723(24) 

1168(2) 

0,1 

0,1 

0,1 

23 55 42.99 

+56 39 14.9 

E 

EB-7 

20.865 

0.560(37) 

0.98(8) 

1.00(13) 

11312(52) 

0,0 

1,0 

1,0 

23 58 27.34 

+56 46 36.4 

E 

EB-8 

18.900 

0.788(9)“ 

0.740(3) 

0.707(3) 

2037(39) 

0,5 

2,3 

5,2 

23 58 29.19 

+56 32 42.7 

RS 

EB-9 

19.981 

0.925(12)“ 

0.701(4) 

0.671(4) 

2647(51) 

1,2 

3,1 

5,2 

23 55 18.41 

+56 43 14.3 

RS 

EB-10 

20.323 

1.025(19) 6 

0.662(8) 

0.630(8) 

2545(79) 

0,0 

1,1 

3,2 

23 57 29.99 

+56 57 34.3 

RS 

EB-11 

20.402 

0.616(24) 

0.89(5) 

0.87(6) 

6784(22) 

0,0 

0,0 

3,2 

23 58 13.37 

+56 45 36.1 

E 

EB-12 

16.724 

0.55(8) 6 

0.94(17) 

0.93(25) 

1589(17) 

0,1 

1,2 

0,0 

23 57 12.86 

+56 31 26.5 

EA 

EB-13 

19.728 

0.562(13)“ 

0.980(23) 

0.988(34) 

6619(59) 

1,1 

3,0 

2,1 

23 55 59.17 

+56 45 14.1 

EA 

EB-14 

17.571 

0.427(10)“ 

1.33(4) 

1.57(7) 

5599(85) 

0,3 

1,2 

3,1 

23 57 18.04 

+56 51 12.0 

EA 

EB-15 

20.619 

1.133(18)“ 

0.620(9) 

0.586(10) 

2378(133) 

1,1 

2,2 

5,1 

23 58 39.38 

+56 36 44.7 

RS or CV 

EB-16 

21.120 

1.853(19) 

0.239(6) 

0.245(4) 

752(4) 

0,1 

1,0 

1,1 

23 56 57.23 

+56 34 03.4 

E 

EB-17 

19.420 

1.298(9)“ 

0.544(5) 

0.507(4) 

1043(14) 

1,0 

2,0 

2,0 

23 56 47.64 

+56 36 28.7 

RS 

EB-18 

17.650 

0.64(9)“ 

0.81(12) 

0.78(14) 

1578(10) 

0,0 

0,1 

0,2 

23 56 01.68 

+56 43 08.3 

RS 

EB-19 

17.368 

0.62(7) 

0.82(10) 

0.79(12) 

1477(2) 

0,1 

0,1 

0,2 

23 56 36.67 

+56 52 43.4 

RS 

EB-20 

20.080 

0.797(18)“ 

0.743(6) 

0.710(5) 

3423(84) 

0,0 

1,0 

2,1 

23 56 11.77 

+56 45 55.6 

RS 

EB-21 

18.984 

1.516(9) 

0.399(6) 

0.379(5) 

526(2) 

0,1 

0,0 

1,0 

23 57 51.16 

+56 42 03.2 

E 

TR-1 

20.703 

0.632(32) 

0.87(5) 

0.84(6) 

7291(30) 

0,0 

0,1 

1,0 

23 57 45.06 

+56 55 36.6 


TR-2 

18.024 

0.467(8) 

1.195(22) 

1.348(37) 

5348(7) 

0,0 

0,0 

1,0 

23 57 24.52 

+56 55 17.6 


TR-3 

19.553 

0.971(9) 

0.679(4) 

0.649(4) 

1995(5) 

1,0 

0,0 

2,1 

23 56 47.11 

+56 51 10.2 



a Colour corrected using the lightcurve model. ^Colour corrected using interpolation of the lightcurve. 


to be constrained for such a system are P, to, i, Rc and a constant 
magnitude mo . Periodic variations in the apparent brightness of the 
star were also accounted for in three different ways, leading to three 
competing planetary transit models: 

h(t) = /o(l - fc{t)) Model 1 

/ 2 (f) = fi(t) ^1 + Asin Model 2 

fs(t) = ~ Ce cos p — j 

' Model 3 

_C h cos(2^1)) 

(16) 

The function /„(/) is the predicted stellar flux at time t for Model 
n, fo is a constant flux value and /c(f) is the fraction of the total 
stellar flux obscured by the companion at time t. The function /i (t) 
is calculated in a numerical fashion by creating a grid for the ob¬ 
served stellar disk and calculating the flux from each grid element 
taking into account the apparent position of the companion at time 
t and the effect of linear limb darkening with u = 0.5. 

Model 1 is therefore appropriate for a star with a constant 
brightness. Model 2 incorporates sinusoidal stellar flux variations 
of semi-amplitude A and phase <f> which do not necessarily have 
the same period as the orbital period of the companion. Such varia¬ 
tions may be present for stars with a lot of star spot activity. Model 
3 incorporates stellar flux variations due to two effects. The first ef¬ 


fect, modelled by the Ce cosine term, is due to the star being tidally 
distorted into an ellipsoidal shape by the companion and rotation- 
ally synchronised. The value of Ce quantifies the semi-amplitude 
of such ellipsoidal flux variations. The second effect, modelled by 
the Cfo cosine term, is due to heating on one side of the compan¬ 
ion caused by irradiation by the star. The value of C^ quantifies the 
semi-amplitude of the heating term. 

If stellar flux variations exist, and they are best modelled by el¬ 
lipsoidal and/or heating terms, then this favours a stellar rather than 
a planetary companion. A planet does not have enough mass to dis¬ 
tort the shape of the star, and neither does a planet emit enough 
radiation to be detectable in the lightcurve. Eclipses with different 
depths also imply a stellar companion since the secondary eclipse 
is caused by occultation of the companion, indicating that the com¬ 
panion is emitting enough radiation to contribute to the observed 
brightness of the system. 

In the following sections, the transit candidates have been 
organised into groups depending on the lightcurve properties, 
and analysed accordingly. Relevant parameters are shown on the 
lightcurve plots. Periods were initially determined using a peri- 
odogram, and then refined during the fitting of the appropriate tran¬ 
sit model. We refer to any transit fit for which all 5 parameters have 
been optimised as a full transit fit. For brevity, the following la¬ 
belling format has been adopted for the plots in these sections: 

<Star No.> - <Plot Code> - <Run(s) To Which The Plot Applies> 

The plot codes are as follows: 

(i) L - Lightcurve 

(ii) PL - Phased lightcurve 

(iii) CPL - Close up of phased lightcurve 
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(a) EB-1 - L & CTFIT - 1999-07 Night 3 (b) EB-2 - CPL & CTFIT - 2000-09 



(c) EB-3 - L & CTFIT - 2000-09 Night 1 



(d) EB-4 - L & CTFIT - 2000-09 Night 8 




(e) EB-5 - L & CTFIT - 2000-09 Night 1 (f) EB-6 - L & CTFIT - 1999-07 Night 6 



(g) EB-7 - L & CTFIT - 2000-09 Night 10 


Figure 6. Eclipsing binaries with undetermined periods. 
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(iv) BL - Binned lightcurve 

(v) CTFIT - Central transit fit 

(vi) FTFIT - Full transit fit 

(vii) CM - Chi squared contour map showing the best fit solution 
with a cross and the 1, 2 and 3 a confidence regions with solid, 
dashed and shorter dashed lines respectively. Annular, grazing and 
no eclipse regions are separated by thick solid lines. 


6.2 Eclipsing Binaries With Undetermined Periods 

In Fig. 6 we present 7 transit candidates for which we were unable 
to determine a period, although we were able to classify them as 
eclipsing binaries. The reason for not being able to determine the 
period was due to either the presence of only one fully/partially 
observed eclipse in the lightcurve and/or cycle ambiguity between 
eclipses. We fitted the best-defined eclipse for each candidate with 
Model 1 keeping the inclination fixed at 90.0° (we call this a central 
transit fit). This fit determines a minimum radius of the companion 
given the eclipse profile, since at lower inclination values the same 
size companion obscures a smaller fraction of the total stellar flux 
due to limb darkening effects and the possibility that the eclipse 
is grazing instead of annular. Flence, at lower inclinations a larger 
companion radius is required to account for the observed eclipse 
depth. 

For each of the 7 transit candidates except EB-3, the central 
transit fit yields a minimum companion radius that is greater than 
0.2Rq. This favours a stellar rather than a planetary companion. 
The lack of out-of-eclipse lightcurve variations leads us to con¬ 
clude that these are eclipsing binaries. For EB-2, a periodogram 
reveals that there are only two possible periods (3.6216±0.0053 d 
and 7.233±0.010 d). We plot the 2000-09 lightcurve folded on the 
shorter period in Fig. 6b. 

In the case of EB-3, the derived minimum companion radius 
of 0.187 ?q is most likely an under estimate since the eclipse is 
possibly deeper than the fit shown in Fig. 6f. Also, the lightcurve 
shows sinusoidal out-of-eclipse variations. Therefore we class this 
system as a RS CVn type eclipsing binary. 


6.3 Eclipsing Binaries Exhibiting Secondary Eclipses 

In Fig. 7 we present 4 transit candidates that exhibit secondary 
eclipses in their lightcurves implying that the companion is lumi¬ 
nous. Fig. 7 shows the folded lightcurve for each transit candidate 
along with the best fit transiting planet model. The lightcurves from 
the different runs are offset vertically (in magnitude) from each 
other in order to highlight any changes in the out-of-eclipse vari¬ 
ations. 


6.4 Eclipsing Binaries Exhibiting Ellipsoidal Variations And 
Heating Effects 

In Fig. 8 we present 3 transit candidates that exhibit ellipsoidal 
variations and heating effects in their lightcurves which immedi¬ 
ately implies that the companion is stellar. Fig. 8 shows the folded 
lightcurve for each transit candidate along with the best fit tran¬ 
siting planet model using Model 3. The lightcurves from different 
runs are offset vertically (in magnitude) from each other in order to 
highlight any changes in the amplitude of the out-of-eclipse varia¬ 
tions. 


6.5 A Possible Long Period Cataclysmic Variable 

EB-15 is a 10.8 h eclipsing binary (Fig. 9) that has round-bottomed 
eclipses lasting 0.1 in phase and orbital modulations that peak near 
phase 0.2. With r' « 20.62 mag and r' — i! « 1.13 mag, the 
star falls close to the cluster main sequence (Fig. 2a). The colour 
index is consistent with a 0.62 Mq K7V star at d = 2.4kpc, a pos¬ 
sible cluster member. Over the three runs the orbital modulations 
increase in amplitude from 0.1 mag to 0.2 mag, while the eclipse 
depth decreases from 0.40 mag to 0.24 mag. 

The orbital modulation could arise from spots on one or both 
stars, though this would require a preferred longitude that remains 
stable over 15 months. The orbital phasing is consistent with that 
of an “orbital hump” that is often seen in quiescent dwarf novae, 
arising from the anisotropic emission of a “hotspot” on the rim of 
an accretion disk where the mass transfer stream from the compan¬ 
ion star feeds material into the disk. The relatively shallow eclipse 
would then imply a moderate inclination so that the donor star 
eclipses only the near rim of the disk and possibly the hotspot. 
The eclipse shape is more symmetric than would be expected for 
eclipses of a hotspot, however, and a hotspot eclipse would become 
deeper rather than shallower as the orbital modulation increased. 
We are unable to decide which interpretation may be correct and 
recommend follow up observations to resolve this ambiguity. In 
any case the eclipse is too deep to be attributed to a planetary tran¬ 
sit. 


6.6 More Eclipsing Binaries 

In this section we present 6 transit candidates exhibiting neither 
easily discernible secondary eclipses nor orbital modulations con¬ 
sistent with ellipsoidal variations and heating effects. Figs. 10 and 
11 show for each star a folded and binned lightcurve with the 
best fit transiting planet model, and an unbinned close-up of the 
folded lightcurve around the primary eclipse along with the best fit 
model. For EB-17 the lightcurves from different runs are offset ver¬ 
tically to highlight changes in the amplitude and phase of the out- 
of-eclipse variations. We rule out the transiting planet model for 
stars EB-16 thru EB-20 because the full transit fit yields a compan¬ 
ion radius greater than 0.2Rq. For EB-21 the full transit fit admits 
Rc < 0.2 Rq but only for periods P < 1.1 d that are ruled out and 
hence the transiting planet model is ruled out also for EB-21. 

EB-16 is an interesting case in its own right. The lightcurve 
shows 4 V-shaped eclipses of depth 0.39 mag and duration 1.3 h 
over the three runs. Figs. 10a and 10b show the lightcurve folded 
on the 1.31 d period. With r' « 21.12 mag and r' — i' ~ 1.85 mag 
we derive a primary 0.24 A/q M5V star at d ~750pc, in front of 
the cluster. The full transit fit reveals that the companion is the same 
size as the primary at 0.2457 ?q and that the eclipses are grazing. 
Hence the period is actually 2.62 d and we classify the system as 
a grazing eclipsing binary consisting of a pair of M5V stars, an 
interesting discovery in that few such systems are known. 

EB-21 is a difficult case in that the lightcurve data show 
one eclipse in the 1999-06 run with only 2 data points during the 
eclipse, and one well sampled eclipse in the 2000-09 run (Fig. 1 If). 
The eclipse is V-shaped of depth 0.2 mag and duration 1.0 h sug¬ 
gesting that it is likely to be a grazing eclipse. With r' ss 18.98 mag 
and r' — i' « 1.52 mag we find that the primary is a 0.40 Mq M2V 
star that lies at only ~530pc. A central transit fit to the 2000-09 
data (Fig. Ilf) yields a minimum companion radius of 0.141 
due to the small size of the primary star (0.38 Rq) from which 
we cannot rule out the transiting planet model. Analysis of the 
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Phase 


(a) EB-8 - PL. BL & FTFIT - 2000-09 (Top), 1999-07 (Middle with 
+0.1 mag offset) & 1999-06 (Bottom with +0.2 mag offset) 



(c) EB-10 - PL, BL & FTFIT - 2000-09 (Top) & 1999-07 (Bottom with 
+0.4 mag offset) 



(b) EB-9 - PL, BL & FTFIT - 2000-09 (Top), 1999-07 (Middle with 
+0.2 mag offset) & 1999-06 (Bottom with +0.4 mag offset) 



Figure 7. Eclipsing binaries exhibiting secondary eclipses. 


shape of the single eclipse in the 2000-09 data is possible due to 
the good time sampling of the observations and such an analysis 
may reveal whether the eclipse is the result of an annular occulta- 
tion by a smaller companion or a grazing occultation by a larger 
companion. Also, we may attempt to predict the orbital period P 
of the planetary companion as a function of the impact parameter 
b = a cos i/R* of the eclipse and subsequently use the lightcurve 
data from the whole run to determine which periods, and hence 
which values of 6, may be ruled out. 

For EB-21 we created a grid for the impact parameter b from 
0.0 to 1.0. For each value of b we fitted Model 1 to the single eclipse 
in the 2000-09 run, using the lightcurve data from the whole run, 
in order to determine a time of mid-transit to, a constant magni¬ 
tude mo, a planetary radius R c and a transit duration At. The chi 
squared Xg C i of the fit was also calculated. The duration of a transit 
event is given by: 


At = 



(17) 


We also have Kepler’s third law: 

2 47r 2 a 3 

P = - 

GM, 


(18) 


We already know M* and R *, and since we have Rc and At as 
functions of b from our fits of the single eclipse, we may use Equa¬ 
tions 17 and 18 to estimate P (or a) as a function of b for the 
transiting planet model. For each value of b we folded the 2000- 
09 lightcurve of EB-21 on the predicted period P using the fit¬ 
ted to, and calculated a new chi squared Xfoid using the fit to the 
single eclipse with period P. In general, if the predicted period is 
such that none of the folded lightcurve data falls during the eclipse, 
then = XecL However, if the predicted period is such that 

some of the folded lightcurve data does fall during the eclipse, then 
Afold XecL ruling out that particular period, impact parameter 
and corresponding eclipse solution. 

Fig. lie shows a plot of Xg C i versus Rc/R* (dashed line), 
which appears constant due to the scale on the y-axis. The con¬ 
tinuous line is a plot of Xfold versus Rc/R*, which clearly shows 
that Xfold Xecl f° r va l ues °f Rc from 0.1417?q (the min- 
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+0.1 mag offset) 



(b) EB-13 - PL, BL & FTFIT - 1999-06, 1999-07 & 2000-09 



+0.1 mag offset) & 1999-06 (Bottom with +0.2 mag offset) 

Figure 8 . Eclipsing binaries exhibiting ellipsoidal variations and heating effects. 


imum companion radius with 6 = 0 and Rc/R* = 0.373) to 
0.227 Rq (b = 0.92 and Rc/R* = 0.600). This is due to the 
fact that the predicted period is less than 1.10 d for these values of 
6. This demonstrates that the transiting planet model is inconsistent 
with our observational data for this star, and hence a stellar com¬ 
panion is favoured. We have classified this system as an eclipsing 
binary. 

6.7 INT-7789-TR-1 

TR-1 exhibits one poorly sampled partially observed eclipse during 
the 1999-07 run and one well sampled fully observed eclipse of 
depth 0.07 mag and duration 6.0 h in the 2000-09 run (Fig. 12b), 
and as a result we were unable to determine a period for the system. 
The star has r' fn 20.70 mag and r' — i' « 0.63 mag from which 
we derive a 0.87Mq late G star primary that lies far behind the 
cluster (d = 7.3kpc). A central transit fit to the 2000-09 data yields 
a minimum companion radius of 0.188-Rq from which we cannot 
rule out the transiting planet model. 

Analysing the single eclipse in the 2000-09 run using the same 


method as for EB-21 in Section 6.6 yields a predicted period of 
7.0 d for b = 0 that increases rapidly with increasing values of b. 
In fact Xfoid = Xecl f° r ^ > 0.09. Fig. 12a shows a plot of 
Xgci versus Rc/R* (dashed line) for 6 = 0 (Rc/R* = 0.223) to 
6 = 1.13 (Rc/R* = 0.500). The minimum value of Xg C i obtained 
is Xecl = 925.0 corresponding to Rc/R* = 0.500, and this is 
marked on Fig. 12a as a horizontal shorter dashed line, along with 
the chi squared values Xg C i = 925.0 +1.0 and Xg C ( = 925.0 + 4.0 
corresponding to the 1 and 2(j confidence levels. 


One can see from Fig. 12a that R c /R* > 0.243 with a lcr 
confidence. This is equivalent to stating that Rc > 0.205-Rq with 
a 1 a confidence. As a result, we can only rule out the transiting 
planet model for this transit candidate at the lcr level and therefore 
further observations are required to confirm the conclusion that this 
system is an eclipsing binary. Fig. 12b shows a plot of the solu¬ 
tion corresponding to the minimum value of Xg C i along with the 
lightcurve data for the night on which the eclipse occurs. This solu¬ 
tion predicts a period of 58+12 d and an inclination of 89.1+0.3°. 
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(a) EB-15 - PL, BL & FTFIT - 2000-09 (Top), 1999-07 (Middle with 
+0.3 mag offset) & 1999-06 (Bottom with +0.6 mag offset) 



+0.3 mag offset) & 1999-06 (Bottom with +0.6 mag offset) 
Figure 9. Possible long period cataclysmic variable. 


Table 7. Star, companion and lightcurve properties for TR-3 as obtained from the various fits detailed in Section 6.9. Column 4 is calibrated r' magnitude. M c 
and Fie are the companion mass and radius respectively. The ratio fc/f* is the flux ratio of the companion to the primary star in the Sloan r' waveband. The 
quantity d is the distance to the system. 


Model 

to - 2451000.0 

P 

m 0 

Am 

At 

i 

M* 

R* 

Me 

Rc 

fc/f* 

d 

x 2 


(d) 

(d) 

(r' mag) 

(mag) 

(h) 

n 

(Mq) 

(Rq) 

(Mq) 

(Rq) 

(Sloan r') 

(pc) 


A 

800.6201(7) 

2.4867(116) 

19.553 

0.086 

1.93 

83.48(6) 

0.701(5) 

0.670(5) 

0.661(30) 

0.629(29) 

0.669 

2788(6) 

725.88 

B 

800.6203(7) 

1.2431(58) 

19.553 

0.070 

1.72 

87.1(1.2) 

0.679(4) 

0.649(4) 

0.000 

0.151(7) 

0.000 

1995(5) 

729.99 

C 

800.6204(7) 

1.2430(58) 

19.553 

0.070 

1.68 

88.8(1.5) 

0.682(4) 

0.651(4) 

0.126(6) 

0.149(6) 

0.0019 

2021(5) 

731.73 


A = Similar Size Stellar Companion. B = Planetary Companion. C = Small Stellar Companion. 


6.8 INT-7789-TR-2 

TR-2 shows a single 0.02 mag eclipse of duration 2.5 h dur¬ 
ing the 2000-09 run (Fig. 12d). With r' « 18.02 mag and 

r' — i 1 « 0.47 mag we find that the primary is a l.lOMpe F star at 
d = 5.3kpc, behind the cluster. A central transit fit to the 2000-09 
data yields a minimum companion radius of 0.174 Rq from which 
we cannot rule out the transiting planet model. 

Applying the same analysis as for EB-21 in Section 6.6 to 
the single eclipse yields a predicted period of 0.62 d for 6 = 0 
that increases slowly with increasing values of b. For b < 0.66, 
Xf 0 id Xecl and f° r b — 0.66, Xfoid oscillates between the 
states Xfoid > Xecl and Xfoid = Xecl- Hence we can be sure 
that b > 0.66, which corresponds to Rc/R* > 0.137. In Fig. 12c 

we plot Xecl versus Rc/R* (dashed line) and Xfoid versu s Rc/R* 
(continuous line) where: 


Xfoid 


Xfoid if b < °- 66 

Xecl if 6 > 0.66 


(19) 


The minimum value of x 2 c j obtained is Xp, = 806.6 correspond¬ 
ing to Rc/R* = 0.137, and this is marked on Fig. 12c as a 
horizontal shorter dashed line, along with the chi squared values 
Xed = 806.6 + 1.0. Xed = 806.6 + 4.0 and xf c \ = 806.6 + 9.0 
corresponding to the 1, 2 and 3<r confidence levels. 

One can see from Fig. 12c that 0.137 < Rc/R* < 0.144 
with a ler confidence. This is equivalent to stating that 
Rc = 0.185!q'qqq Rq. Flence the conclusion at the lcr level 


is that this is a possible transiting planet in orbit around a 1.20Mg 
F star that merits follow up observations. Fig. 12d shows a plot 
of the solution corresponding to the minimum value of Xg C i along 
with the lightcurve data for the night on which the eclipse occurs. 
This solution predicts a period of 1.8±1.3 d and an inclination of 
81.9±1.6°. 


6.9 INT-7789-TR-3 

TR-3 exhibits two fully observed eclipses and one partially ob¬ 
served eclipse during the 2000-09 run. The eclipses have a depth of 
0.07 mag and duration 1.7 h with a period of 1.24 d (Figs. 12e, g 
& h). The star has r' « 19.55 mag and r' — i 1 « 0.97 mag from 
which we derive a 0.68 A/q primary star of spectral type K5V that 
lies slightly in front of the cluster (d = 2.0kpc). A full transit fit 
to the 2000-09 lightcurve data yields a best fit companion radius of 
0.151±0.007i?Q consistent with the radius of a transiting planet 
(Fig. 12f). This solution is reported in Table 7 under Model B, and 
the °f the fit is 729.99. 

Flowever, there are two other models for TR-3 that should be 
considered to see if they produce a better f 2 value for the fit to the 
lightcurve data. It is possible that the companion is a smaller and 
less luminous star than the primary star that produces a secondary 
eclipse which is not visible in the lightcurve folded at the ~ 1.24 day 
period, and it is also possible that the companion is a star of similar 
size and luminosity to the primary star and that the system actually 
has an orbital period of ~2.48 days. In order to test these models 
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(a) EB-16 - PL, BL & FTFIT - 1999-06, 1999-07 & 2000-09 
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(b) EB-16 - CPL & FTFIT - 1999-06, 1999-07 & 2000-09 
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(c) EB-17 - PL, BL & FTFIT - 2000-09 (Top), 1999-07 (Middle with 
+0.2 mag offset) & 1999-06 (Bottom with +0.4 mag offset) 
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(d) EB-17 - CPL & FTFIT - 2000-09 (Top), 1999-07 (Middle with 
+0.2 mag offset) & 1999-06 (Bottom with +0.4 mag offset) 
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Figure 10. EB-16, EB-17 and EB-18. 
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(a) EB-19 - PL, BL & FTFIT - 2000-09 






0.45 0.5 

Planet To Star Radius Ratio 

. ,2 


(e) EB-21 - Xq C \ (dashed line) & Xf 0 id (continuous line) versus R c /R* 
2000-09 



Figure 11. EB-19, EB-20 and EB-21. 
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(a) TR-1 - Xed (dashed line) versus Rc/R* - 2000-09 
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(b) TR-1 - L & Single Eclipse Fit - 2000-09 Night 9 



(c) TR-2 - Xq C \ (dashed line) versus R c / R * & Xfoid ( con_ 
tinuous line) versus Rc/R* - 2000-09 
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(g) TR-3 - PL, BL & FTFIT - 2000-09 
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(h) TR-3 - CPL & FTFIT - 2000-09 
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Secondary Star To Primary Star Radius Ratio Secondary Star To Primary Star Radius Ratio 

(a) TR-3 - CM - 2000-09 (b) TR-3 - CM - 2000-09 



(c) TR-3 - PL, BL & Eclipsing Binary Fit - 2000-09 (d) TR-3 - PL, BL & Eclipsing Binary Fit - 2000-09 
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(e) TR-3 - CPL & Eclipsing Binary Fit - 2000-09 
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(f) TR-3 - CPL & Eclipsing Binary Fit - 2000-09 



Figure 13. Eclipsing binary fits for TR-3. 
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for TR-3 we have developed an eclipsing binary model based on 
the same assumptions as for the star and planet system presented 
in Section 6.1 except that we assume that the companion is now 
luminous and massive, and that our theoretical main sequence rela¬ 
tionships adopted in Section 4.3 apply to the companion. 

The eclipsing binary model has five parameters to optimise: 
orbital period P, time of mid-eclipse to, orbital inclination i, com¬ 
panion to primary star radius ratio R c /R* and a constant magni¬ 
tude mo. We fitted this model to the lightcurve of TR-3 by calcu¬ 
lating the \ 2 for a grid in i and R c / R* . For each value of R c /R „, 
we had to recalculate the distance d to the system, the values of 
M* and R , for the primary star, and the mass of the companion 
M c . This was done by constructing a theoretical binary main se¬ 
quence for the current value of Rc / R, and then finding the dis¬ 
tance d such that this model passes through the position of TR-3 
in the colour-magnitude domain. The initial value of P was either 
1.24 d or 2.48 d corresponding to the smaller or similar size stel¬ 
lar companion models respectively. Table 7 reports the results of 
these fits. Fig. 13 shows a chi squared contour map, a folded and 
binned lightcurve with the best fit eclipsing binary model, an un¬ 
binned close-up of the folded lightcurve around the primary eclipse 
along with the best fit model and another unbinned close-up around 
the secondary eclipse along with the best fit model. The left hand 
column of diagrams in Fig. 13 applies to the case of the smaller 
stellar companion and the right hand column of diagrams in Fig. 13 
applies to the case of the similar size stellar companion. 

Table 7 shows that the best model for TR-3 is the eclipsing bi¬ 
nary model with a similar size stellar companion since this model 
attains the smallest \ 2 of 725.88. All three models require exactly 
5 parameters to be optimised and hence we calculate a likelihood 
ratio of ~7.8 for the eclipsing binary model with a similar size 
stellar companion compared to the transiting planet model, and we 
calculate a likelihood ratio of ~18.6 for the eclipsing binary model 
with a similar size stellar companion compared to the eclipsing bi¬ 
nary model with a smaller stellar companion. Finally, we calcu¬ 
late a likelihood ratio of ~2.4 for the transiting planet model com¬ 
pared to the eclipsing binary model with a smaller stellar compan¬ 
ion. Flence our conclusion is that this system is most likely to be 
a grazing eclipsing binary with period 2.49 d consisting of a K4V 
star primary and a K5V star secondary that lies at d = 2.8kpc, 
slightly behind the cluster. Flowever, further observations will be 
required confirm this conclusion and categorically rule out the tran¬ 
siting planet model. 


7 CONCLUSIONS 

In the search for our transit candidates we have developed an accu¬ 
rate, efficient and fast photometry pipeline employing the technique 
of difference image analysis. Raw data from the telescope are pro¬ 
cessed by the pipeline with minimal user input in order to directly 
produce lightcurves and colour magnitude diagrams. This is espe¬ 
cially important considering the high quantity of data that may arise 
from a transit survey. 

Our analysis of the colour magnitude diagrams by including 
the treatment of extinction for the open cluster NGC 7789 has al¬ 
lowed us to assign a model-dependent mass, radius and distance to 
each star. Such information is vital in the subsequent analysis of 
the transit candidates since it allows a direct estimate of the com¬ 
panion radius. We detected 24 transit candidates which warranted 
a detailed analysis of their lightcurves and we were able to deter¬ 
mine periods for 14 of these candidates. Of the 10 candidates with¬ 


out periods, we could rule out the transiting planet model for 7 of 
them by determining the minimum companion radius and for an¬ 
other one by predicting the orbital period. For INT-7789-TR-1, it 
was only at the ler level that we could rule out the transiting planet 
model based on the shape of the best eclipse. For INT-7789-TR- 
2 we found a companion radius of 0.1852qqqq-Rq (1.812q'qq7?j) 
based on the analysis of the best eclipse. Follow-up observations 
(see below) will be required for both of these candidates in order 
confirm that INT-7789-TR-1 is an eclipsing binary and in order to 
determine the nature of INT-7789-TR-2. 

For the 14 transit candidates with well determined periods, we 
could rule out the transiting planet model for 4 of them from the 
detection of previously disguised secondary eclipses, and for 3 of 
them from the observation that the out-of-eclipse lightcurve data 
exhibit ellipsoidal variations and heating effects. One of the can¬ 
didates is possibly a new cataclysmic variable with a long period 
(10.8 h) which could be a cluster member, worthy of follow up ob¬ 
servations in its own right. All of the 8 above mentioned candidates 
plus another 5 may be ruled out as having planetary companions by 
considering that the companion radius obtained from the full transit 
fit is greater than 0.2 Rq . For INT-7789-TR-3, none of the above ar¬ 
guments may be used to rule out the transiting planet model. How¬ 
ever, on application of an eclipsing binary model to the lightcurve 
we find that the model consisting a pair of grazing K dwarf stars 
is ~7.8 times more likely than the transiting planet model. This 
is by no means a definitive conclusion that INT-7789-TR-3 is an 
eclipsing binary since there is a non-negligible probability that the 
transiting planet model is still valid. Follow-up observations will 
be required to confirm that INT-7789-TR-3 is the type of eclipsing 
binary that we predict in this paper. 

Future photometric observations of the three transit candidates 
for which we could not rule out the transiting planet model with 
confidence should consist of time series observations in two differ¬ 
ent filters. Eclipsing binary status may be confirmed by the observa¬ 
tion of different eclipse depths in different filters since a planetary 
transit is an achromatic event. INT-7789-TR-1 and INT-7789-TR-2 
also require the observations of multiple eclipses in order to deter¬ 
mine their period and whether they exhibit secondary eclipses or 
not. If these follow up photometric observations still allow the pos¬ 
sibility that the transiting planet model is valid, then radial velocity 
observations may be used to place an upper limit on the mass of 
the orbiting companion, hopefully low enough to rule out a stel¬ 
lar or brown dwarf companion. The fact that these candidates are 
so faint ( r' « 20.7 mag for INT-7789-TR-1, r' « 18.0 mag for 
INT-7789-TR-2 and r' « 19.6 mag for INT-7789-TR-3) makes it 
very unlikely that radial velocity observations with 10-m class tele¬ 
scopes will achieve the accuracy required to determine the actual 
mass of the companion iCharbonneail2003l ). 

From a simple signal-to-noise argument presented in Section 
4.5, we expected to detect ~2 transiting hot Jupiters. At most we 
have detected 3 transiting hot Jupiters, but our analysis of these 
candidates shows that this is very unlikely. Follow up observations 
will most likely show that our candidates are eclipsing binaries, 
which means that our transit survey will have produced a null re¬ 
sult. We are currently modelling in more detail the number of hot 
Jupiters that we expected to detect as a function of planetary radius, 
orbital period, star mass and detection threshold via Monte Carlo 
simulations. The results of these simulations are consistent with our 
simple estimate of the expected number of detections, and by com¬ 
bining this information with the results of further observations on 
the three remaining transit candidates, we will be able to estimate 
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the hot Jupiter fraction of the cluster field, a very important result 
for the testing of star and planetary formation theories. 
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